Attractive Fermi gases with unequal spin populations in highly elongated traps 
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We investigate two-component attractive Fermi gases with imbalanced spin populations in trapped 
one dimensional configurations. The ground state properties are determined within local density 
approximation, starting from the exact Bethe-ansatz equations for the homogeneous case. We 
predict that the atoms are distributed according to a two-shell structure: a partially polarized 
phase in the center of the trap and either a fully paired or a fully polarized phase in the wings. The 
partially polarized core is expected to be a superfluid of the FFLO type. The size of the cloud as 
well as the critical spin polarization needed to suppress the fully paired shell, are calculated as a 
function of the coupling strength. 
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Recent experimental studies of polarized Fermi gases 
near a Feshbach resonance have attracted a great deal 
of interest [H0]- Since a finite spin polarization desta- 
bilizes the Bardeen-Cooper-Schrieffer (BCS) s-wave su- 
perconductivity, these strongly interacting systems are 
indeed promising candidates to search for exotic pairing 
mechanisms. 

So far, the experiments have mainly explored three- 
dimensional (3D) trapped configurations. The MIT data 
[H show that in the presence of a shallow trap the po- 
larized Fermi gas phase-separates into three concentric 
shells: a fully paired superfluid core, an intermediate par- 
tially polarized normal shell and, finally, a fully polarized 
phase at the edge of the cloud. The Rice experiments 
have instead been performed on more anisotropic cigar- 
shaped configurations and the obtained measurements 
suggest a somewhat different scenario Q. A plausible 
explanation is that in the Rice experiments the applied 
trapping potential in the radial directions is tight enough 
that the gas cannot be considered as locally homogeneous 
and qualitatively new effects are seen which cannot be 
accounted for by the 3D local density approximation. 

In this Letter, we study the limiting case where the ra- 
dial confinement is so strong that only the lowest trans- 
verse mode is populated and the gas is kinematically one- 
dimensional (ID). We show that the ground state of the 
ID polarized Fermi gas in a shallow axial trap exhibits a 
two-shell structure, in sharp contrast with the 3D case. 
Our analysis is based on the exact equation of state for a 
ID homogeneous Fermi gas with arbitrary strong attrac- 
tive interactions and the ID local density approximation. 

Strongly interacting ID Fermi gases have already been 
investigated experimentally using ultracold atoms in op- 
tical lattices [3(. The possibility to create spin imbal- 
anced configurations opens new prospects to detect the 
Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) state [i] [j| un- 
der clean experimental conditions. In particular, it was 
recently argued by Yang [6j that the ground state of ID 
homogeneous attractive Fermi gases with unequal spin 
populations is the ID analogue of the FFLO state. For 



trapped configurations, this implies that the partially po- 
larized phase of the ID gas is a superfluid of the FFLO 
type and the experimentally measurable cloud sizes and 
density profiles provide indirect signatures of this unusual 
state. This is clearly exciting, because FFLO supercon- 
ductivity has been of interest to condensed matter physi- 
cists for decades and, more recently, it has attracted the 
attention of the particle physics community [3] ■ 

Consider a two-component atomic Fermi gas in the 
presence of an axially symmetric harmonic potential 
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where m is the atom mass and r± = y ' x 2 + y 2 . We 
require the anisotropy parameter A = lu z /lu±_ to be suf- 
ficiently small that only the lowest transverse mode is 
populated. At zero temperature, this condition implies 
that the Fermi energy €f] associated with the longitudi- 
nal motion of the majority | component is much smaller 
than the distance between energy levels in the transverse 
direction, epf <g; hai±, or, equivalently, A <C 1/iVf . Un- 
der this assumption, Bergeman et al. [8[ have shown that 
the scattering properties of atoms are well described by 
an effective ID contact interaction giuS(x), with 
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where a 3 o is the 3D scattering length, a± = y / h/muj± 
is the transverse oscillator length and A = 1.0326. 
The effective ID interaction is attractive, g\u < 0, if 
l/ a 3-D < A/a±. By tuning the 3D scattering length 
near a Feshbach resonance, the coupling constant ((2|) 
can be changed adiabatically from the weakly interact- 
ing regime, g±o —* — , to the strong coupling limit, 

giD — > -00. 

The effective ID Hamiltonian Hid acting on the 
fermions is given by H 1D = H® D + X)i=i Vho(Zi), where 
Vho = to(jj 2 z 2 /2 is the longitudinal trapping potential 
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with iV = iVf +iVj_ being the total number of atoms. The 
Hamiltonian ([3]) is exactly soluble by Bethe's ansatz 0. 

Unpolarized ID Fermi gases, in both homogeneous and 
trapped configurations, have been investigated in Ref.flpj 
and Ref. ll|, as an exactly solvable model for the BEC- 
BCS crossover. The homogeneous case with a finite spin 
polarization has been recently studied in Ref. [12| . In the 
following we investigate the magnetic properties of spin 
polarized ID attractive Fermi gases in harmonic traps at 
zero temperature. We first focus on the homogeneous 
case by setting Vho = 0. The effects of the axial confine- 
ment will be discussed later, using the ID local density 
approximation (LDA). 

For fixed values of the linear number densities n-\ = 
iVf /L and ri\ = NJL, where L is the size of the system, 
the ground state energy E of H° D is given by [9(] 



and AE = f\?it 2 a-iD[n\{2n^ - nj + 8nj(n T - njj 3 ]/12m. 
The above result is in agreement with Ref. [l2j ■ Equation 
(|6|) shows that in the strongly interacting regime the sys- 
tem is a mixture of a Tonks gas of 2Vj_ diatomic molecules 
of mass 2m and binding energy es = h 2 /ma\ D , and an 
ideal Fermi gas formed by the unpaired f fermions. 

For a fixed value of the total density n = n.f + n^, 
the density difference s = n-j- — can vary in the range 
< s < n. For s = 0, the ground state of H® D is fully 



paired with a gap in the spin excitation spectrum 11 1 
whereas in the opposite limit s — n, the system is a fully 
polarized gas of f fermions. For < s < n, the gas is 
partially polarized and it is expected to be a superfluid 
of the FFLO type @. 

From Eq.|4]), we calculate the chemical potential 
p = d(E/L)/dn and the effective magnetic field h = 
d(E I L)/ds, which are related to the chemical poten- 
tials p a of the two components by p = (/if + P{)/2 and 
h = - Ml)/ 2 - From Eq.® it follows that 



h = e B F h (n,s), 
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The spectral functions a(X),p(k) in Eq.((4|) are solutions 
of the coupled integral equations 
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where the kernels are given by -Ki(A, A') = l/[ir + tt(X' — 
A) 2 ] and K 2 (X,k) = 2/[tt + 47r(fc - A) 2 ]. Here B and 
Q are non-negative numbers related to the number den- 
sities by n^diD = 2 J „ a(X)dX and n^aiD — nj^aic = 

2 p{k)dk, respectively, and a±n = —2h 2 /mgiD is the 
effective ID scattering length. The coupling strength is 
controlled by the parameter naio, where n = n-\ + n± is 
the total density: the weakly interacting case corresponds 
to nan) 3> 1, whereas the strongly attractive regime is 
achieved for na\D -C 1. 

The leading terms in the weak coupling expansion for 
the energy can be calculated by perturbation theory, 
yielding E/L — h 2 7r 2 (n^ + n^)/6m — 2h 2 n^ni/(maio)- 
In the strongly interacting regime, corresponding to small 
values of B and Q, we expand the spectral functions <r(A) 
and p(k) around A = and k = 0, respectively, retaining 
up to quadratic terms. The coefficients of the expansion 
can be easily obtained from Eq.©. From Eq.((4]), we then 
obtain E = E° + AE, where the leading term is given by 
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where and Fh are universal functions of the dimen- 
sionless number densities h = na\u and s = sauj. The 
energetic stability of the mixture in the partially polar- 
ized phase allows us to use p, and h as new coordinates, 
since the Jacobian of the transformation ((7| is positive 
0. 

By setting s = in Eq. ([7]) we obtain the critical curve 
Pc/es shown in FigQ] and corresponding to the phase 
boundary with the fully paired phase. We see that p, c di- 
verges logarithmically for vanishing magnetic fields and 
decreases as h increases, reaching its minimum value 
(x c = —€b/2 for h — > es/2~. The corresponding asymp- 
totic behavior p c ~ — h is plotted with a dashed line in 
the inset of Fig[Tl where we have magnified the region 
of the phase diagram near the point 0=(es/2, — es/2). 
Similarly to what happens for the attractive Hubbard 
model [lj], the critical magnetic field h c = h(n,0) is re- 
lated to the energy gap A in the spin sector for the unpo- 
larized case by h c = A/2. Since A measures the strength 
of interactions, h c increases by decreasing the density 
n (or, equivalently, the chemical potential p,) reaching 
h c = (b/2 in the strongly attractive limit p — > —eb/2, 
as shown in FigfTJ 

By setting s = n in Eq. |7J) we obtain the second critical 
curve p s /eB plotted in FigQ] and corresponding to the 
phase boundary with the fully polarized phase. In this 
particular case, the saturation field h(n, s = n) and the 
corresponding chemical potential p s = p(n, s = n) can be 
calculated in closed form from Eqs ([4]) and ([5]), yielding 

h = 2e B [Q 2 (l-arctan(2Q)/7r)+(2Q-arctan(2Q)+7r)/47r] 

(8) 

and p s = —h + 2Q 2 £B, where Q = 7mai£>/2. In the weak 
coupling limit, corresponding to h eb, the chemical 
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FIG. 1: Ground state magnetic properties for the homo- 
geneous system. The dotted vertical line corresponds to 
h — es/2. The arrow shows the trajectory of the local chem- 
ical potential in a trap. Inset: zoom of the phase diagram 
near the point 0=(es/2, — es/2): the asymptotic behaviors 
are shown with dashed lines. 



potential diverges as fj, s = h — Allies + (4/7r 2 — l/2)e_g 
whereas in the strong coupling regime it approaches O as 
fx a = -e B /2 + (2 v / 2/37r % /eB)(/i - e B /2) 3 / 2 . This latter 
behavior is also shown in the inset of FigQ] with dashed 
line. 

In FigHJ the fully paired and the fully polarized phases 
are limited from below by a shaded area representing the 
vacuum state, i.e. the state of zero total density. For 
h < eb/2, the minimum value of the chemical potential 
[X = (A v is given by /i„ = — es/2, whereas for h > es/2, we 
find fi v — —h, since fj, + h = fi-\ is the chemical potential 
of the majority component. 

We see from Fig[T] that for a fixed value of the mag- 
netic field, only two phases are allowed for the ID system: 
partially polarized and fully paired if h < es/2 or par- 
tially polarized and fully polarized if h > es/2. This is 
clearly different from the 3D case, where fi c is an increas- 
ing function of the magnetic field 15[ and thus all three 
phases are allowed for a given h. 

The effects of a shallow harmonic potential Vho(z) can 
be taken into account via LDA, which is applicable pro- 
vided the size of the cloud is large compared to the axial 
oscillator length a z = y^h/mtOz, implying N 3> 1. Under 



this assumption, the density profiles n a (r) for the two 
components a =t,l arc obtained by imposing the lo- 
cal equilibrium condition \i a [n-f (z) , + Vho(z) = /v> 
where A*o- [^f? n l\ are the corresponding chemical poten- 
tials of the homogeneous system and are constants 
fixed by the normalization N a = J n a (z)dz. This is 



vi n l( z ), n l( z )} = M° -V ho (z), 



(9) 



where fi°,h° = (fi® ± /z°)/2. Equation (JOj) shows that 
the local mean chemical potential decreases by increasing 
the distance from the center of the trap while the iocai 
magnetic field remains constant throughout. In the phase 
diagram of Fig[TJ this corresponds to follow the vertical 
line h = h° downwards starting from the point fi = fi , 
as shown by the arrow. 

We immediately see that trapped gases with spin im- 
balanced components have a two-shell structure, with an 
inner partially polarized phase in the center of the trap. 
The outer shell is either fully paired, if h° < es/2, or 
fully polarized, if h° > es/2. Hence, the fully paired 
and the fully polarized phases are mutually exclusive in 
a trap. For the special case h° = £_b/2, shown with dot- 
ted line in Fig[TJ the partially polarized phase extends to 
the edge of the cloud. 

In order to relate fjP, h° to the total number N = iV-p + 
N± of atoms and the spin polarization P = (JVj —N±)/N, 
we need to calculate the density profiles, starting from 
Eqs ([3]) and ©. The normalization conditions can be 
expressed in terms of the dimensionless chemical poten- 
tial fl° = h°/cb and the dimensionless densities n and s, 
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where p, v = fi v /eB- Equation (|10[) shows that, in the 
presence of the trap, the coupling strength is given by 
Na\ D /a 2 z , which is large in the weakly interacting regime 
and vanishes in the strongly attractive limit. 

The condition h° = es/2 yields the critical spin po- 
larization P = P c needed to suppress the fully paired 
shell; for P > P c the outer shell is always fully po- 
larized. In Fig[2] we plot the critical spin polarization 
as a function of the parameter Na\ D /a 2 z (solid line). 
We see that P c increases going towards the strong cou- 
pling regime where it saturates to P c = 0.2. This value 
is surprisingly small compared to the critical polariza- 
tion P 3D - 0.77 found [lj in 3D Fermi gases at reso- 
nance (a^D — > oo). The corresponding density profiles 
in the strongly attractive limit, where Eq.© holds, are 
given by n x (z) = (8m/^ 2 ) 1/2 (/ + e B /2 - V ho {z)f/\ 
with fi° + es/2 = Huj z Ni/2, and nj-(z) = ni(z) + 
(2m/ 7 r 2 fi 2 ) 1 /2( M o + h o _ T//io ( z ))i/2 j with + h o _ 

(TVj — Ni)hhj z . The inclusion of the next order correc- 
tion AE yields P c = 1/5(1 - 0.364 v'iVa^/a 2 ). This 
asymptotic behavior is shown in Fig[2]with dashed line. 

For P / P C) the cloud exhibits a two-shell structure, 
with the radius Ri n of the inner shell fixed by the con- 
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FIG. 2: Critical spin polarization as a function of the interac- 
tion parameter Na\ D /a 2 z . The behavior in the strong coupling 
regime is shown with dashed line. The arrow corresponds to 
the asymptotic limit P c — 0.2. 
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FIG. 3: Size of the cloud (left) and radius of the inner shell 
(right) as a function of the spin polarization for different val- 
ues of the parameter Naf D /al — oo, 10, 1, 0.1, 0. The non 
monotonic behavior is a signature of the two-shell structure. 



dition fj, - Vho(Rin) = l*c(h°), for h° < e B /2, and 
fjP — Vh (Ri n ) = n s (h ), otherwise. The size of the 
cloud, where the total density vanishes, is instead given 
by the condition /j, — Vho(R^) = Hv{h )- Both quantities 
are plotted in Fig[3] as a function of the spin polariza- 
tion P for different values of the parameter Na 2 D /a 2 . 
In the absence of interaction (top curve), we find i?| = 
s/1 + Pa 7 N lj2 and R in = Vl^PazN 1 / 2 , showing that 
R-\{Rin) increases (decreases) as P increases. For finite 
interactions, we see that R-\ decreases as P increases 
until P = P c , and it increases for larger values of the 
spin polarization. In particular, in the strong coupling 
regime, we find Rt = ^(1 - P)/2a z N 1 / 2 for P < 1/5, 
and i? T = VZPcizN 1 / 2 for P > 1/5. This non monotonic 
behavior is a signature of the outer shell. Indeed, an in- 
crease in the spin polarization implies that the radius of 
the density profile for the minority component decreases, 
as Ni decreases. For P < P c , the profiles for the two 
components match in the fully paired outer shell, and 
thus also i?| decreases. 

The inner radius (right panel) shows the opposite be- 
havior: Ri„ = for zero spin polarization and it increases 
rapidly until P — P c , where Ri n — i?| and the fully 
paired shell disappears. For higher values of the polar- 
ization, Ri n decreases and finally vanishes at P — 1. In 
particular, in the strong coupling regime, we find Ri n = 
VTPazN 1 / 2 , for P < 1/5, and R m = ^(1 - P)/2a z N 1 / 2 
for P > 1/5. 

The above scenario for a ID polarized Fermi gas can be 
investigated in current experiments with ultracold atoms 
under a tight radial confinement, e.g. 40 K samples in an 
array of ID tubes created with optical lattices [3j. In 



order to reach a strongly degenerate regime, the ther- 
mal energy k B T should be small compared to the Fermi 
energy, implying k B T <C NHuj z . This condition can be 
satisfied by taking the typical experimental parameters 
Q uj z /2tt ~ 250 Hz for the axial frequency, mean num- 
ber of particles (per tube) N - 100 and T - 50nK. The 
spin imbalance can be easily induced by a radio- frequency 
sweep and the density profiles measured by absorption 
imaging techniques. 

In conclusion, we have presented an exact solution of 
an experimentally accessible ID model of a polarized 
Fermi gas at arbitrary strong coupling. We have found 
that, differently from the 3D case, the trapped gas phase 
separates in a two-shell structure, with a FFLO super- 
fluid core surrounded by either a fully paired or a fully 
polarized phase, depending on the value of the spin polar- 
ization. Experimental signatures of this novel and exotic 
state are a reduced value of the critical spin polarization 
P c and a non monotonic behavior of the size of the cloud 
as a function of the spin polarization. 
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